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Abstract 

We use QCD sum-rules to predict ground state masses for pseudoscalar ( J PC = (T + ) charmonium and bottomonium 
hybrids. We find that the inclusion of a six-dimensional gluon condensate contribution is needed to stabilize the 
analyses. For the charmonium hybrid, we find a mass of (3.82 ± 0.13) GeV; for the bottomonium hybrid, we find a 
mass of (10.64 + 0.19) GeV. We comment on possible phenomenological implications concerning the F(3940). 
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Over the past decade, more than a dozen new 
charmonium-like states, the XYZ resonances, have been 
discovered (see yj] for a recent review). Few of these 
states can be easily accommodated with a conventional 
charmonium meson interpretation [2]. There are dis- 
crepancies between observations and calculations re- 
garding masses and widths; certain J PC sectors seem 
to be overpopulated — the vector (1 ) states, in partic- 
ular; and electrically charged resonances have perhaps 
been seen. Naturally, there has been considerable spec- 
ulation that some of these new states may lie outside of 
the constituent quark model. 

Motivated by these findings, we use QCD sum-rules 
to predict ground state masses of both pseudoscalar 
(0~ + ) charmonium and bottomonium hybrids. The first 
application of sum-rules to heavy hybrids was done 
in ||2r0|. Therein, a variety of J PC quantum numbers 
were investigated; however, in some cases (such as 1" 
and 0~ + ), the sum-rules were unstable and the result- 
ing mass estimates were deemed unreliable. Recently, 
some of these sum-rules have been re-analysed 16J48J]. 
In the vector analysis of [6], it was seen that the in- 
clusion of a six-dimensional gluon condensate contribu- 
tion stabilized the sum-rules and led to a reliable ground 
state mass prediction. Therefore, in an effort to update 
the original pseudoscalar work in a similar fashion, we 
also include a six-dimensional gluon condensate contri- 
bution. 
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As in J3], we define a current 



(i) 



where Q represents a heavy quark (charm or bottom) 
operator, and we define a corresponding two-point func- 
tion 

Tl, v (q) = if d\ «**<0| Tj,{x)j v {Q) |0) (2) 



-gJn {V) (q 2 ) + ^n (S) (q 2 ). (3) 



(The S stands for scalar; the V for vector.) Then, the 
longitudinal projection of 11^ 



q- 



(4) 



serves as a probe of heavy, pseudoscalar hybrids. 

We compute II(S) using the operator product expan- 
sion (OPE) with contributions from perturbation theory, 
the four-dimensional gluon condensate 



(aG 2 ) = {aG%G a n 

= (7.5 + 2.0) x 10~ 2 GeV 4 

as well as the six-dimensional gluon condensate 



~ 3 C 3 > = {g i f ahc G%G bv p G c P> 1 ) 

(8.2 ± 1.0)GeVV(? 2 >, 



(5) 



(6) 



and we write 

W) - n ^ 2 ) + + n W)- (7) 

September 20, 2012 



The numerical values in (O and (|6]l are extracted from 
heavy quark systems [9]. At leading-order in a, the 
relevant Feynman diagrams are shown in Figure 
Divergences are handled using dimensional regulariza- 




cut discontinuity in the complex z-plane extending from 
Z = 1 to z — > 00 ; as such, it is readily seen from (TST> — (TTOb 
that fl^ D has a branch cut extending from q 2 = 4m 2 to 
q 2 — > 00 as expected. 

Due to its analytic structure and asymptotic be- 
haviour, I1(S) satisfies the following dispersion relation: 




Figure 1: Feynman diagrams included in the OPE of II(S). 

tion in D = 4 + 2e dimensions followed by MS- 
renormalization. The Wilson coefficients of both gluon 
condensate terms are calculated using fixed-point gauge 
methods (see 01 1 |, | l 2ll. for example). We use the pro- 
gram TARCER 11311 to apply the two-loop integral re- 
currence relations of 11141 Il5ll which significantly re- 
duces the number of integral formulae needed to ob- 
tain exact results. The required two-loop integrals are 
in lfl6ll ; the required one-loop integrals are in Jl7ll . Ex- 
panding in e and omitting polynomials in q 2 as they ul- 
timately provide no contribution to the sum-rules, we 
find 



^ 2 ) = ^-\to {z - 1) ^- 21z + 10 ) 



x 3*2(1,1,1; 3/2, 3; z) 
+ ^z(8z 3 + 8z 2 +29z-10) 



X 3*2(1,1, 2;5/2,4;z) 



. (8) 



n (S)V)= lk z(2z+ 1 )2*'i(l. 1; 5/2; z)(«G 2 >, (9) 



and 



W) = 



1 



(2z 2 - 2z + 1) 



384tt 2 (z- l) 2 

x 2 *i(l,l;5/2;z) 

+ (10z 2 - 20z + 7) 



? 3 G 3 > (10) 



where z = and where p F q {- ■■ ; z) are general- 
ized hypergeometric functions (see [18], for example). 
For q = p-l, the function p F q (- ■ ■ ; z) has a branch 
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4m 2 



+ e 2 n[ S) (0) + -e 4 n' ( 's)(°) + ze 6n w(0) (ID 



where Q 2 = -q 2 . To eliminate unwanted polynomials 
in q 2 and to suppress the high-energy (continuum) con- 
tributions from Imll(S) to the integral on the right-hand 
side of (fTTT i. we apply the Borel transform 



~k#~ T(N) \dQ 2 

N/Q 2 =t 



(12) 



to obtain 



'All Feynman diagrams are drawn using JaxoDraw| 10]. 



1 k 1 f°° 

-s [( - e 2 ) n (S) (e 2 )] = - rV T imn (S) (o dt. 

(13) 

On the left-hand side of ( TT3l . we use (TTb — (TTOb to approx- 
imate ri(S); on the right-hand side, we employ a single 
narrow resonance plus continuum model 

Imn (S )(0 = f 2 6 (t -M 2 ) + 6(t- s ) ImlT^ D (?) (14) 

where sq is the continuum threshold parameter and / 
represents a hadronic coupling. Moving the threshold 
contribution from the right- to the left-hand side of ( TT3l , 
we arrive at the Laplace sum-rules B19L 12011 

Z k (r,s )^[(-Q 2 ) k ^(Q 2 )} 

_I r t k e-' T lmn^ D (t)dt (15) 

n Js 

=f 2 M 2k e- M2 \ 

By exploiting a well-known relationship between the 
Borel transform and the inverse Laplace transform 112 111 , 
we can simplify the remaining Borel transformed term 
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in ( [TBI ) (see B22I1 . for example). Doing so gives 



For bottom quarks 



Am 1 r Io/4m , r / 

£ ( T , So ) = ^- J ^ 4 '" ™ [imn^f (4m 2 x) 

+IiriIT G 3 G (4m 2 x)] dx 

Am 2 r So/4m2 

— e - 4m ""Imn^ G (4m 2 x)^ 
4m 2 (g 3 G 3 > _ 4m ,/ 



+ lim 

;7-»o + 



and 



128tt 2 V^7 



Xi (t, s ) = -7-X0 (t, so) 
or 



(16) 



(17) 



where, for z > 1 , 



imn-y ) = 



Vz - 1 Vz 



120tt 2 z 2 

x (30 - 1 15z + 166z 2 + 8z 3 + 16z 4 ) 
-15(-2 + 9z-16z 2 + 16z 3 )log(Vz :r T+ Vi) 



nn/ ,\ m 2 Vz — 1 1 

Imnf^ 2 ) = — (2z+ l)^L^<aG 2 ), 



(18) 



(19) 



and 



imnf s f( ? 2 ) 



1 vfh - 

256ttz(z - l) 2 ^ 

x (2z 2 - 2z + l)<g 3 G 3 >. (20) 



We note that the portion of dTOb that has a double pole 
at z = 1 ultimately contributes to both ( TToT ) and ( fT7b ; 
hence, Imll^ 1 alone is insufficient to formulate the 
sum-rules. A full calculation of IT^ 00 is required. 

In ([Tol l and ( TTTl i. the parameters m and a are MS- 
scheme running quantities evaluated at a scale yiz. For 
charm quarks, 



a(jx) 



a (M T ) 



1 + 



25a(M r ) 
12k 



—,a(M T ) = 033 (21) 



_ / a(jj) 
m c Qi) = m c 

\a(m c ) 



12/25 



Th c = (1.28 + 0.02) GeV. 

(22) 



a(jj) 



a(M z ) 



1 + 



23q-(M z ) 



a(M z ) = 0.118 (23) 



Wi(jU) = OTf, 

\ar(m 6 ) 



T2^ lo g(ft) 

12/23 

, m 6 = (4.17 + 0.02) GeV. 



(24) 



In (I2"T1>- (|2"4"1 >. the quark mass parameters are from 
23 425ll . the t and Z masses are from j26Tl . and 
a(M T ) and a(Mz) are from B27I1 . Furthermore, 
renormalization-group improvement sets = -^= lEsTl . 
From ( fT5b , it follows that 



£l(T, So) = m2 

' £o(t, s ) 



(25) 



The range of acceptable r-values is determined by de- 
manding that, on the left-hand side of d25l ). the con- 
tinuum contributes less than 30% and the condensates 
contribute less than 10%. Then, so and M are cho- 
sen to optimize (l25l l over this range. See [3] for fur- 
ther details. In the charmonium hybrid case, we find 
s = 23 GeV 2 and M = (3.82 ± 0.13) GeV; in the 
bottomonium hybrid case, we find s = 140 GeV 2 and 
M = (10.64 + 0.19) GeV. The uncertainties in these 
mass predictions are dominated by the uncertainties in 
the quark mass parameters (l22l i and (124-b and the un- 
certainties in the six-dimensional gluon condensate (0. 
Plots of the left-hand side of d2"5b for various values of 
so are shown in Figures [2] and [3] for charmonium and 
bottomonium hybrids respectively. Inclusion of a six- 
dimensional gluon condensate contribution does indeed 
seem to stabilize the sum-rule analysis. 

VX1/X0 (GeV) 
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Figure 2: The ratio -£^ CD (t, so) /X? CD (t, so) for hybrid charmonium 
is shown as a function of the Borel scale 1 /t for the optimized value 
so = 23 GeV - (solid curve). For comparison the ratio is also shown 
for so = 28 GeV 2 (upper dotted curve) and so = 19 GeV 2 (lower 
dotted curve). The uppermost dashed curve represents the so + * °° 
limit corresponding to the bound M < 3.96 GeV. Central values of 
the QCD parameters have been used. 
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Figure 3: The ratio X^ CD (t, sq) /jS CD (t, sq) for hybrid bottomo- 
nium is shown as a function of the Borel scale 1/r for the optimized 
value sg = 140 GeV 2 (solid curve). For comparison the ratio is also 
shown for so = 155 GeV 2 (upper dotted curve) and sq = 116 GeV 2 
(lower dotted curve). The uppermost dashed curve represents the 
sq — > co limit corresponding to the bound M < 10.84 GeV. Central 
values of the QCD parameters have been used. 



The 7(3940) (likely the same resonance as the 
Z(3915)) seen by both the Belle JH and BaBar Q 
collaborations is a C = + isosinglet with a mass of 
3.915 GeV. Unfortunately, at present, its parity is un- 
known and its spin is uncertain: either J = or J — 
2. As such, the 7(3940) could be a pseudoscalar, but 
additional work is needed to establish its J p assign- 
ment. In i29ll . this particle was already touted as a 
hybrid candidate based on both its production mecha- 
nism and its observed decay mode. The 7(3940) is seen 



in B decays which are thought to be particularly well- 
suited to hybrid production D3lt- Also, the T(3940) has 

not been seen to decay through kinematically allowed 

— (*) 

DD modes, an observation consistent with a flux 
tube model-inspired selection rule which suppresses hy- 
brid decays to pairs of S -wave mesons Il32l434ll . If the 
T(3940) is ultimately shown to have J pc = 0~ + , then 
our mass prediction of (3.82 ±0.13) GeV would provide 
additional support in favour of a charmonium hybrid in- 
terpretation. 
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